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We study thermoelectric devices where a single 18-annulene molecule is connected to metallic 
zigzag graphene nanoribbons (ZGNR) via highly transparent contacts that allow for injection of 
evanescent wave functions from ZGNRs into the molecular ring. Their overlap generates a peak in 
the electronic transmission, while ZGNRs additionally suppress hole-like contributions to the ther- 
mopower. Thus optimized thermopower, together with suppression of phonon transport through 
ZGNR-molecule-ZGNR structure, yield the thermoelectric figure of merit ZT ~ 0.5 at room tem- 
perature and 0.5 < ZT < 2.5 below liquid nitrogen temperature. Using the nonequilibrium Green 
function formalism combined with density functional theory, recently extended to multiterminal 
devices, we show how the transmission resonance can also be manipulated by the voltage applied to 
a third ZGNR electrode, acting as the top gate covering molecular ring, to tune the value of ZT. 
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Thermoelectrics transform temperature gradients into 
electric voltage and vice versa. Although a plethora of 
thermoelectric-based energy conversion and cooling ap- 
phcations has been envisioned, their usage is presently 
limited by their small efficiency [T]. Careful tradeoffs 
are needed to optimize the dimensionless figure of merit 
ZT — S^GT/k quantifying the maximum efficiency of a 
thermoelectric conversion because ZT contains unfavor- 
able combination of the thermopower S, average temper- 
ature T, electrical conductance G and the total thermal 
conductance k = KcI + ^ph (including contributions from 
electrons Hci and phonons Kph)- The devices with ZT > 1 
are regarded as good thermoelectrics, but ZT > 3 is re- 
quired to compete with conventional generators jl] . 

The major experimental efforts to increase ZT have 
been focused on suppressing the phonon conductivity us- 
ing either complex (through disorder in the unit cell) 
bulk materials [5] or bulk nanostructured materials [3]. 
A complementary approach engineers electronic density 
of states to obtain a sharp singularity P] near the Fermi 
energy which can enhance the power factor S^G. 

The very recent experiments |1] and theoretical stud- 
ies ISHS] have ignited exploration of devices where a single 
molecule is attached to metallic |6j or semiconducting f7\ 
electrodes, so that dimensionality reduction and possi- 
ble strong electronic correlations [5] allow to increase 
S concurrently with diminishing Kph while keeping the 
nanodevice disorder-free |10j . For example, creation of 
sharp transmission resonances near the Fermi energy Ep 
by tuning the chemical properties of the molecule and 
molecule-electrode contact can substantially enhance the 
thermopower S which depends on the derivative of the 
conductance near Ep. At the same time, the presence 
of a molecule in the electrode-molecule-electrode het- 
erojunction severely disrupts phonon propagation when 
compared to homogenous clean electrode. 




FIG. 1: (Color online) Schematic view of the pro- 
posed ZGNR| 18-annulene|ZGNR three-terminal heterojunc- 
tion. The contact between the source and drain 8-ZGNR 
(consisting of 8 zigzag chains) metallic electrodes and a ring- 
shaped molecule is made via 5-membered rings of carbon 
atoms (dark blue), while the electrodes are attached to atoms 
1 and 10 of the molecule. The third electrode is coupled as an 
"air-bridge" top gate, made of ZGNR as well, covering only 
molecular ring at the distance 5.3 A. The two-terminal ver- 
sion of the device assumes that such top gate is absent. The 
hydrogen atoms (light yellow) are included to passivate the 
edge carbon atoms. 

In this Letter, we exploit a transparent contact be- 
tween metallic zigzag graphene nanoribbon (ZGNR) elec- 
trodes and a ring-shaped 18-annulene molecule to pro- 
pose two-terminal (i.e., top gate absent in Fig. [T]) and 
three-terminal devices illustrated in Fig. [l] Their ther- 
moelectric properties are analyzed using the very recently 
developed nonequilibrium Green function combined with 
density functional theory formalism for multiterminal 
nanostructures (MT-NEGF-DFT) [TTl US]- The high 
contact transparency allows evanescent modes from the 
two ZGNR electrodes to tunnel into the molecular re- 
gion and meet in the middle of it (when the molecule is 
short enough [13]), which is a counterpart of the well- 
known metal induced gap states in metal-semiconductor 
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FIG. 2: (Color online) Physical quantities determining ther- 
moelectric properties of the two-terminal device shown in 
Fig. [l] (without the top gate electrode): (a) zero-bias 
electronic transmission %i{E); (b) the corresponding ther- 
mopower S at two different temperatures; (c) thermoelectric 
figure of merit ZT vs. energy at two different temperatures; 
and (d) ZT vs. temperature at three different energies. 



FIG. 3: (Color online) (a) The phonon transmission function 
Tpi^{uj) and (b) the corresponding thermal conductance Kph 
for an infinite 8-ZGNR and the two-terminal device shown 
in Fig. [T] (assuming absence of the third top gate electrode) 
whose source and drain electrodes are made of semi-infinite 
8-ZGNR. 



Schottky junctions. This effect can induce a large peak 
(i.e., a resonance) in the electronic transmission function 
near Ef [Fig. [2];a)], despite the HOMO-LUMO energy 
gap of the isolated molecule. The enhancement of the 
thermopower [Fig. [sjb)] due to transmission resonance 
around E > Ep and suppression of hole-like contribu- 
tion (i.e., negligible transmission around E < Ep) to 
S, together with several times smaller phonon thermal 
conductance (Fig. [s]) when compared to infinite ZGNR, 
yields the maximum room-temperature ZT ^ 0.5 in the 
two-terminal device [Fig. [2][c),(d)]. Furthermore, we dis- 
cuss how a third top gate ZGNR electrode covering the 
molecule, while being separated by an air gap in Fig. [T] 
can tune the properties of the transmission resonance via 
the applied gate voltage thereby making possible further 
enhancement of ZT (Fig.[4|. 

Among the recent theoretical studies of molecular ther- 
moelectric devices via the NEGF-DFT framework [BJ [7] , 
most have been focused |B] on computing the ther- 
mopower S, with only a few [7] utilizing DFT to obtain 
forces on displaced atoms and then compute Kph- More- 
over, due the lack of NEGF-DFT algorithms for multiter- 
minal nanostructures, the possibility to tune thermoelec- 
tric properties of single-molecule devices via the usage of 
the third electrode has remained largely an unexplored 
realm 0. We note that the recent proposal |2| for the 
two-terminal molecular thermoelectric devices with so- 
phisticated combination of a local chemical tuning of the 
molecular states and usage of semiconducting electrodes 
has predicted much smaller ZT ~ 0.1 at room tempera- 



ture. In addition, our 0.5 < ZT < 2.5 at E-Ep = -0.02 
eV (which can be set by the backgate electrode cov- 
ering the whole device [23]) in the temperature range 
T = 30-77 K is much larger than the value achieved in 
conventional low-temperature bulk thermoelectric mate- 
rials g]. 

The recent fabrication of GNRs with ultrasmooth 
edges [H], where those with zigzag edges are insulating 
at very low temperatures due to edge magnetic ordering 
which is nevertheless easily destroyed above > 10 K [15j . 
has opened new avenues for highly controllable molecular 
junctions with a well-defined molecule-electrode contact 
characterized by high transparency, strong directionality 
and reproducibility. This is due to the fact that strong 
molecule-GNR tt-tt coupling makes possible formation of 
a continuous 7r-bonded network across GNR and orbitals 
of conjugated organic molecules [13]. Unlike the metallic 
carbon nanotubes (CNTs) employed experimentally [IB] 
to generate such networks [13j, GNRs have planar struc- 
ture appropriate for aligning and patterning. 

The early experiments [TB] on CNT| molecule |CNT het- 
ero junctions have measured surprisingly small conduc- 
tances for a variety of sandwiched molecules. The first- 
principles analysis of different setups reveals that this is 
due to significant twisting forces when molecule is con- 
nected to CNT via, e.g., 6-membered rings [13] ■ There- 
fore, to keep nearly parallel and in-plane configuration 
(hydrogen atoms of 18-annulene slightly deviate from the 
molecular plane) of our ZGNR|18-annulene| ZGNR junc- 
tion, we use a 5-membered ring [IF in Fig. [l] to chemi- 
cally bond ZGNR to annulene. 
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The carbon atoms of a ring-shaped 18-annulene 
molecule can be connected to ZGNR electrodes in con- 
figurations whose Feynman paths for electrons traveling 
around the ring generate either constructive or destruc- 
tive quantum interference effects imprinted on the con- 
ductance [T7]. For example, a 7r-electron at Ep entering 
the molecule in setup (1,10) shown in Fig.[T]has the wave- 
length kp/2d {d is the spacing between carbon atoms 
within the molecule), so that for the two simplest Feyn- 
man paths of length 9d (upper half of the ring) and 9c? 
(lower half of the ring) the phase difference is 0. Note that 
the destructive quantum interference [T7| would form an 
additional dip (i.e., antiresonance) within the main trans- 
mission peak in Fig. [2][a) generated by the injection of 
overlapped evanescent modes [12]. The effect of antires- 
onance on the thermopower S for gold|18-annulene|gold 
junctions has been studied in Ref. [T5] as a possible sen- 
sitive tool to confirm the effects of quantum coherence 
on transport through single-molecule junctions. 

The computation of quantities entering ZT for real- 
istic single-molecule junctions requires quantum trans- 
port methods combined with first-principles input about 
atomistic and electronic structure to capture charge 
transfer in equilibrium (which is indispensable to obtain 
correct zero-bias transmission of, e.g., carbon-hydrogen 
systems [12]), geometrically-optimized atomic positions 
of the molecular bridge including molecule-electrode sep- 
aration in equilibrium, and forces on atoms when they 
are perturbed out of equilibrium. The state-of-the- 
art approach that can capture these effects, as long 
as the coupling between the molecule and the elec- 
trodes is strong enough to ensure transparent contact 
and diminish Coulomb blockade effects [H dO] , is NEGF- 

DFT [laiin]. 

The technical details of the construction of the 
nonequilibrium density matrix via NEGF-DFT for mul- 
titerminal devices are discussed in Ref. [TT]. Our 
MT-NEGF-DFT code utilizes uhrasoft pseudopoten- 
tials and Perdew-Burke-Ernzerhof (PBE) exchange- 
correlation functional. The localized basis set for DFT 
calculations is constructed from atom-centered orbitals 
(six per C atom and four per H atom) that are opti- 
mized variationally for the electrodes and the central 
molecule separately while their electronic structure is ob- 
tained concurrently. 

In the coherent transport regime, the NEGF post- 
processing of the result of the DFT self-consistent loop 
expresses the zero-bias electron transmission function be- 
tween the left (L) and the right (R) electrodes as: 

%m = Tr {Tl{E)G{E)Tr{E)G\E)] . (I) 

The matrices Tl^r{E) — i[T,L,R{E) — j^{E] account 
for the level broadening due to the coupling to the elec- 
trodes, where 11l^ji{E) are the self-energies introduced 
by the ZGNR electrodes [T^]. The retarded Green func- 
tion matrix is given by G = [i?S — H — — S^]"^, 



where in the local orbital basis {(pi} Hamiltonian matrix 
H is composed of elements Hij = ((/)i|i?Ks|'^j) (^KS is 
the effective Kohn-Scham Hamiltonian obtained from the 
DFT self-consistent loop) and the overlap matrix S has 
elements Sij = {4>i\4>j). 

The transmission function Eq. ([I]) obtained within the 
NEGF-DFT framework allows us to compute the follow- 
ing integrals |21] 

C30 

K^{^.)^l J dE%,{E){E~f^r(^~^^^y (2) 

— OO 

where f{E,fi) = {l + cxp[{E-fi)/kBT]}~'^ is the 
Fermi-Dirac distribution function at the chemi- 
cal potential /i. The knowledge of i^„(/i) finally 
yields all electronic quantities in the expression for 
ZT: G = e^Ko[y)] S = Ki{y)/[eTKo{ti)]\ and 
K,, = {K2{^i) - [K^{tx)YlKo{ti)}/T. 

The phonon thermal conductance is obtained from the 
phonon transmission function T-ph{Lo) using the corre- 
sponding Landauer-type formula for the scattering re- 
gion (molecule + portion of electrodes) attached to two 
semi-infinite electrodes: 

oo 

'^P'^ = 2^fc^ j ^"""^P''(" \eWfe^^-l)^ - 



The phonon transmission function 7^h(w) can be calcu- 
lated from the same Eq. ([!]) with substitution H — >• K 
and S w^M, where K is the force constant matrix 
and M is a diagonal matrix with the atomic masses. We 
obtain the force constant matrix using GPAW, which is 
a real space electronic structure code based on the pro- 
jector augmented wave method [12] • The electronic wave 
functions are expanded in atomic orbitals with a single- 
zeta polarized basis set, and PBE exchange-correlation 
functional is used. The whole scattering region, which in- 
cludes 27 layers of ZGNR electrodes, is first relaxed to a 
maximum force of 0.01 eV/Aper atom. Subsequently, we 
displace each atom / by Qia in direction a = {a;, y, z} to 
get the forces Fia,,jfs on atom J ^ I in direction /3. The 
elements of K are then computed from finite differences, 
K/o,j/3 = [FjfiiQic) - Fjfi{-Qi^)]/2Qi^. The intra- 
atomic elements are calculated by imposing momentum 
conservation, such that Kjajjs ~ —Yij^jKiajp- 

Figure [2][a) shows the zero-bias electronic transmission 
Tc\{E) for the two-terminal version of the device in Fig.jl] 
where the peak near Ep is conspicuous. Additionally, 
the suppression of the hole- like transmission \Tc\{E) — ?> 
around E < Ep\ evades unfavorable compensation [7] 
of hole-like and electron-like contributions to the ther- 
mopower. This is due to the symmetry of the valence 
band propagating transverse mode in GNR semi-infinite 
electrode which changes sign at the two carbon atoms 
closest to the molecule for the geometry (molecule con- 
nected to the middle of the GNR edge) in Fig. [I] We 
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FIG. 4: (Color online) Physical quantities determining ther- 
moelectric properties of the three-terminal device shown in 
Fig. [l] as a function of the applied gate voltage Vg: (a) zero- 
bias electronic transmission Tci[E, Vg); (c) the corresponding 
thermopower S at two different temperatures; (c) thermoelec- 
tric figure of merit ZT vs. temperature aX Vg = Q; and (d) 
ZT vs. Vg at two different temperatures. 



emphasize that these features of Tc\{E) are quite insen- 
sitive to the details of (short) conjugated molecules, and 
since they are governed by the ZGNR Bloch states, they 
are impervious to the usual poor estimates of the band 
gap size and molecular energy level position in DFT. 

The maximum value of the corresponding ther- 
mopower S plotted in Fig. [2jb) is slightly away from 
E — Ep = and it is an order of magnitude larger 
than the one measured on large-area graphene [23] or 
in molecular junctions with gold electrodes More- 
over, the interruption of the infinite ZGNR by a molecule 
acts unfavorably to phonon transmission 7ph(w), thereby 
generating three times smaller Kph at room temperature 
when compared in Fig. [3] to the thermal conductance of 
an infinite 8-ZGNR. Figures [2|^c),(d) demonstrate that 
the interplay of large S'^G and reduced Kph for the pro- 
posed two-terminal device yields the room-temperature 
ZT - 0.5 around E- Ef = 0. 

The introduction of the narrow air-bridge top gate 
electrode in Fig. [l] which is positioned at the distance 
5.3 A (ensuring negligible tunneling leakage current into 
such third ZGNR electrode) away from the two-terminal 
device underneath while covering only the molecular 
ring, makes possible tuning of the transmission resonance 
shown in Fig. |4|^a). Even in the absence of any applied 
gate voltage (Vg = 0), ZT vs. temperature plotted in 
Fig. is notably modified when compared to the cor- 
responding functions in Fig. [2jd) for the two-terminal de- 
vice. This stems from slight hybridization of the top gate 
and molecular states. The narrowing of the transmission 



peak around Ep due to the application of negative gate 
voltage enhances the thermopower in Fig. |4][b), thereby 
increasing ZT above its value at Vg — 0, which can be 
substantial at low temperatures as shown in Fig. |4|^d) . 

In conclusion, we predict that a single conjugated 
molecule attached to metallic GNR electrodes, where 
the transparent molecule-GNR contact allows evanes- 
cent modes to penetrate from the electrodes into the 
HOMO-LUMO molecular gap generating a transmission 
resonance, can act as an efficient thermoelectric device. 
Our first-principles simulations suggest that its figure 
of merit can reach ZT ^ 0.5 at room temperature or 
0.5 < ZT < 2.5 below liquid nitrogen temperature, which 
is much higher than ZT found in other recent proposals 
for molecular thermoelectric devices J7|. Moreover, in- 
troduction of the third air-bridge top gate covering the 
molecule can change the sharpness of the transmission 
resonance via the application of the gate voltage, thereby 
opening a path toward further optimization of ZT. We 
anticipate that much higher ZT could be achieved by 
testing different types of molecules to reduce Kph since the 
power factor S^G is already optimized in our device by 
the usage of GNR electrodes which generate molecular- 
level-independent transmission resonance while the sym- 
metry of propagating modes in GNRs lifts the compen- 
sation of hole-like and electron-like contributions to S. 
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